Doping dependence of chemical potential and entropy in hole- and electron-doped 

high-T c cuprates 
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We examine the thermodynamic properties of the hole- and electron-doped cuprates by using the 
t-t'-t"-J model. We find that the chemical potential shows different doping dependence between the 
hole and electron dopings. Recent experimental data of the chemical potential shift are reproduced 
except for lightly underdoped region in the hole doping where stripe and/or charge inhomogeneity are 
expected to be important. The entropy is also calculated as a function of the carrier concentration. 
It is found that the entropy of the electron-doped system is smaller than that of the hole-doped 
systems. This is related to a strong antiferromagnetic short-range correlation that survives in the 
electron-doped system. 
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High-T c superconductivity emerges with carrier dop- 
ing into insulating cuprates. The carrier is either an 
electron or a hole. The phase diagram exhibits an 
asymmetry between the electron and hole dopings: in 
the electron-doped cuprate Nd2- a; Ce a ;Cu04 (NCCO), an 
antiferromagnetic (AF) phase remains up to the con- 
centration x ~ 0.15, while in the hole-doped cuprate 
La 2 _2;Sr2Cu04 (LSCO) the AF phase disappears with 
an extremely small amount of x. 1 Remarkable differ- 
ences of the electronic states between the two materi- 
als have been observed in several experiments. Inelas- 
tic neutron scattering experiments showed the presence 
of incommensurate spin structures in LSCO but not in 
NCCOi^ The optical conductivity exhibits a gaplike fea- 
ture at around 0.2 eV in the AF phase of NCCO^ but 
there is no such feature in LSCO. From angle-resolved 
photoemission experiments, it is clearly observed that 
hole carriers doped into the parent Mott insulators first 
enter into the (±7r/2,±7r/2) points in the Brillouin zone, 
but electron carriers are accommodated at (±7r,0) and 
(0,±7r)i» The doping dependence of the core-level pho- 
toemission also shows different behaviors of the chemical 
potential shift between NCCO and LSCO£ It is inter- 
esting that even in hole-doped cuprates the chemical po- 
tential shift strongly depends on materials^: the shift is 
larger in Bi 2 Sr 2 Cai_ a; (Pr,Er) 2; Cu 2 8+ . y (BSCCO) than 
in LSCO. This indicates the difference of the electronic 
states among the hole-doped cuprates. 

In previous studies^ we showed that the t-J model 
with long-range hoppings t' and t" explains the differ- 
ences of the inelastic neutron scattering, optical conduc- 
tivity, and angle-resolved photoemission data between 
hole- and electron-doped cuprates. In this paper, we ex- 
amine the thermodynamic properties of the hole- and 
electron-doped cuprates by using the same model. A 
finite-temperature version of the Lanczos method for 
small clusters is employed to calculate the thermody- 
namic properties under the grand canonical ensemble. 
The calculated chemical potential shows a different de- 
pendence on the carrier concentration between the hole 
and electron dopings. The experimental data-*& are 



nicely reproduced except for a lightly underdoped re- 
gion in the hole doping where stripe and/or charge inho- 
mogeneities are expected to play important roles. The 
entropy is also calculated as a function of the carrier 
concentration. It is found that the entropy in NCCO 
is smaller than the entropies of LSCO and BSCCO. A 
strong AF short-range correlation that survives in the 
electron-doped system is the origin of the small entropy. 

The t-J Hamiltonian with long-range hoppings, termed 
the t-t'-t"-J model, is 
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where the summations (i, j) lst , 
over first-, second- and third- nearest-neighbor pairs, re- 
spectively. No double occupancy is allowed, and the rest 
of the notation is standard. In the model, the difference 
between hole and electron carriers is expressed as the sign 
difference of the hopping parameters^ t > 0, if < 0, and 
t" > for hole doping, and t < 0, if > 0, and t" < for 
electron doping. The ratios t'/t and t" /t are taken to be 
material dependent: (t'/t, t"/t)=(-0.M, 0.23) for NCCO 
and BSCCO, and (-0.12, 0.08) for LSCO£ J/\t\ is taken 
to be 0.4. In order to examine the thermodynamic prop- 
erties of the model, we use a finite-temperature version of 
the Lanczos method 9 for a square lattice with \fl8x\fl8 
sites under periodic boundary conditions^ The chemi- 
cal potential /i and the entropy density s are calculated 
under the grand-canonical ensemble. 

Figure 1 shows the carrier concentration x for the t-t'- 
t"-J model with different parameter values as a function 
of n at T = J/4 = 0.1 \t\. The data for the t-J model 
are consistent with previous reports With increasing 
magnitudes of t' and t" in the hole carrier side (from t-J 
to LSCO and BSCCO), the slope of the ivs/i curves at 
a small concentration region x < 0.2 becomes weak. The 
derivative of x with respect to // is proportional to the 
charge compressibility k oc —dx/d/i. The variation of the 
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FIG. 1: Carrier concentration a; vs chemical potential /i for 
several parameter sets of the t-t'-t"-J model with 18 sites. 
T I \t\ = 0.1 and J / \t\ = 0.4. The hopping parameters are 
(t,t',*")=(l,0,0) for t-J (dotted line), (1,-0.12,0.08) for 
LSCO (solid line), (1,-0.34,0.23) for BSCCO (dashed line), 
and (-1,0.34, -0.23) for NCCO (dot-dashed line). 



slope thus means that k decreases with increasing t' and 
t", i.e., in the order t-J, LSCO, and BSCCO. The fact 
that the charge fluctuation weakens with the increase of 
long-range hoppings is consistent with the tendency that 
the phase-separated region in the x vs J/t phase diagram 
at zero temperature is reduced with increasing i'^i On 
the other hand, the slope in NCCO is similar to that of 
t-J, indicating that the charge fluctuation is as strong 
as that in t-J. Comparing NCCO and BSCCO, both of 
which have the same magnitude of \t'\ and \t"\, we can 
see a remarkable difference in the doping dependence of 
fi between the electron- and hole-doped systems. 

In Fig. 2, the experimental data of the chemical poten- 
tial shift Afi (Refs. 5 and 6) are replotted for the sake of a 
comparison with our theoretical results. The experimen- 
tal shift is measured from ijl at the lowest concentration 
in each panel. In the lightly doped region, almost no 
change of A/j, is observed in the hole-doped materials, 
while in NCCO A/i is proportional to x. It is also inter- 
esting in the experimental data that the two hole-doped 
materials LSCO and BSCCO exhibit different behaviors: 
the change of A/i at around x ~ 0.1 is larger in BSCCO 
than in LSCO. 

In order to compare our results in Fig. 1 with the 
experimental A^, we use the following procedures: (i) 
We assume that \t\ is material independent and has 
|i| = 0.35 eV£ (ii) We examine a carrier concentra- 
tion x* where /x is almost temperature independent, 9 
and fit fi at x* to an experimentally expected A/x 
at x*: (a;*,A/i)=(0.18,-0.1 eV), (0.24,-0.4 eV), and 
(0.16,0.2 eV) for LSCO, BSCCO, and NCCO, respec- 
tively. We find in Fig. 2 that our results nicely re- 
produce the global features of the experimental A/j, in 
both electron- and hole-doped materials. In particular, 
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FIG. 2: Chemical potential shift Afi vs carrier concentra- 
tion x for both hole- and electron-doped systems. The solid, 
dashed, and dot-dashed lines represent the calculated results 
for the parameters of LSCO, BSCCO, and NCCO, respec- 
tively. The value of \t\ is set to be 0.35 eV for all systems. 
J/\t\ = 0.4 and T/\t\ = 0.1. The circles and triangles de- 
note the experimental data of LSCO and NCCO taken from 
Ref. 5, respectively, and the squares are the BSCCO data 
from Ref. 6. The experimental shift is measured from fi at 
the lowest concentration examined in the experiments. 



the calculated data show a good agreement with the ex- 
perimental data at around x ~0.1 in both LSCO and 
BSCCO. Therefore, the different chemical potential shifts 
between the two materials can be attributed to the dif- 
ference in the long-rang hoppings. Although the global 
agreement with experiment is satisfactory, we find re- 
markable deviations from the experimental data in the 
lightly underdoped regions of the hole-doped systems 
(x < 0.1), where experimentally d/j,/dx ~ (k — > oo). 
One of the possible origins of the deviations might be the 
difference of the temperature between the measurements 
(T - 80 K) and the calculations (T = 0.1 \t\ = 350 mcV 
~ 400 K). If we were able to reduce T by an increase 
in the system size, the deviations would become small 
because of the enhancement of However, it is un- 
likely that k diverges independently of the magnitude 
of t' and t" . In this context, we hope for experiments 
at higher temperatures. Another possibility that can ac- 
count for the deviations might be the presences of stripes 
and/or charge inhomogeneity that are experimentally re- 
portediiSii^ Since we have no evidence of the charge in- 
homogeneity in our calculations, it may be necessary to 
calculate /x in the presence of external potentials that 
induce such an inhomogeneity. This would be a future 
problem to be solved. 

Figure 3 shows the doping dependence of the entropy 
density s at various temperatures for the three param- 
eter sets of the t-t'-t"-J model. At Tj \t\ = 0.4, there 
is no remarkable difference among LSCO, BSCCO, and 
NCCO, particularly in the underdoped region. With 
decreasing temperature from T / \t\ — 0.4 to 0.1, s in 
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FIG. 3: Entropy density s vs carrier concentration x for 
both hole- and electron-doped systems with 18 sites. The 
solid, dashed, and dot-dashed lines represent the results for 
the parameters of LSCO, BSCCO, and NCCO, respectively. 
J I \t\ = 0.4. The three sets of the three lines correspond 
to the data at T j |f |=0.1, 0.2, and 0.4 from the bottom to 
the top. The circles shows experimental results for LSCO at 
T = 320 K, taken from Ref. 14. 



NCCO is strongly suppressed as compared with LSCO 
and BSCCO. We find by examining the temperature de- 
pendence of the spin correlation^ that the suppression of 
s in the underdoped region is correlated with the devel- 
opment of the AF short-range order. This is easily under- 
stood because the AF order reduces the entropy coming 
from the spin degree of freedom. On the other hand, the 
difference of s between hole- and electron-doped systems 
in the overdoped region is probably due to the difference 
of the density of states. It is desirable to confirm the 
suppression of s in electron-doped systems experimen- 
tally. In Fig. 3, we also plot s measured for LSCO at 



T = 320 Ki4 The agreement with the calculated LSCO 
data at Tj \t\ = 0.1 is qualitatively good, but not quan- 
titatively satisfactory. There are several reasons for the 
disagreement: (i) uncertainties in the conversion of a the- 
oretical T into a realistic one and in the experimental de- 
termination of s, (ii) a finite-size effect in our calculations 
that is seen as the dip of s at x ~ 0.22, which comes from 
relatively large sparseness of the low-energy levels in the 
18-site four-hole system, and (iii) more plausibly the ef- 
fect of the stripe and/or charge inhomogeneity discussed 
above, by which s in the underdoped region is expected 
to be reduced^ 

In summary, we have examined the thermodynamic 
properties of the hole- and electron-doped cuprates by 
using the t-t'-t"-J model. The calculated chemical po- 
tential shows different behaviors between the hole and 
electron dopings. The experimental data of the chemical 
potential shift are explained by taking into account the 
material dependences of t' and t", except for the lightly 
underdoped region in the hole doping where the stripe 
and/or charge inhomogeneities are expected to be im- 
portant. The entropy is also calculated as a function of 
the carrier concentration. It is found that the entropy 
of the electron-doped system is smaller than that of the 
hole-doped ones. This is related to a strong AF short- 
range correlation that survives in the electron-doped sys- 
tem. To confirm this, specific heat measurements in the 
electron-doped materials are desired. 
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